Visualization of vascular injuries in extremity trauma.

A tandem of particle-based computational methods is adapted to simulate injury and hemorrhage in the human body. In order to ensure anatomical fidelity, a three-dimensional model of a targeted portion of the human body is reconstructed from a dense sequence of CT scans of an anonymized patient. Skin, bone and muscular tissue are distinguished in the imaging data and assigned with their respective material properties. An injury geometry is then generated by simulating the mechanics of a ballistic projectile passing through the anatomical model with the material point method. From the injured vascular segments identified in the resulting geometry, smoothed particle hydrodynamics (SPH) is employed to simulate bleeding, based on inflow boundary conditions obtained from a network model of the systemic arterial tree. Computational blood particles interact with the stationary particles representing impermeable bone and skin and permeable muscular tissue through the Brinkman equations for porous media. The SPH results are rendered in post-processing for improved visual fidelity. The overall simulation strategy is demonstrated on an injury scenario in the lower leg.

Current Methods and Advances in Simulation of Hemorrhage after Trauma.

As animal models fall out of favor, there is demand for simulators to train medical personnel in the management of trauma and hemorrhage. Realism is essential to the development of simulators for training in the management of trauma and hemorrhage, but is difficult to achieve because it is difficult to create models that accurately represent bleeding organs. We present a simulation platform that uses real-time mathematical modeling of hemodynamics after hemorrhage and trauma and visually represents the injury described by the model. Using patient-specific imaging, 3D-mesh representations of the liver were created and merged with an anatomically accurate vascular tree. By using anatomically accurate representations of the vasculature, we were able to model the cardiovascular response to hemorrhage in a specific artery. The incorporation of autonomic tone allowed for the calculation of bleeding rate and aortic pressures. The 3D-mesh representation of the liver allowed us to simulate blood flow from the liver after trauma. For the first time, we have successfully incorporated tissue modeling and fluid dynamics with a model of the cardiovascular system to create a simulator. These simulations may aid in the creation of realistic virtual environments for training.

GPU-Based Parallelized Solver for Large Scale Vascular Blood Flow Modeling and Simulations.

Cardio-vascular blood flow simulations are essential in understanding the blood flow behavior during normal and disease conditions. To date, such blood flow simulations have only been done at a macro scale level due to computational limitations. In this paper, we present a GPU based large scale solver that enables modeling the flow even in the smallest arteries. A mechanical equivalent of the circuit based flow modeling system is first developed to employ the GPU computing framework. Numerical studies were employed using a set of 10 million connected vascular elements. Run-time flow analysis were performed to simulate vascular blockages, as well as arterial cut-off. Our results showed that we can achieve ~100 FPS using a GTX 680m and ~40 FPS using a Tegra K1 computing platform.

Optimization Integrator for Large Time Steps.

Practical time steps in today's state-of-the-art simulators typically rely on Newton's method to solve large systems of nonlinear equations. In practice, this works well for small time steps but is unreliable at large time steps at or near the frame rate, particularly for difficult or stiff simulations. We show that recasting backward Euler as a minimization problem allows Newton's method to be stabilized by standard optimization techniques with some novel improvements of our own. The resulting solver is capable of solving even the toughest simulations at the [Formula: see text] frame rate and beyond. We show how simple collisions can be incorporated directly into the solver through constrained minimization without sacrificing efficiency. We also present novel penalty collision formulations for self collisions and collisions against scripted bodies designed for the unique demands of this solver. Finally, we show that these techniques improve the behavior of Material Point Method (MPM) simulations by recasting it as an optimization problem.

Cardiovascular blood flow analysis under normal and open injury conditions.

The aim of this paper is to enable a simulation tool for cardiovascular blood flow in order to better understand normal and hemorrhage conditions. A second order partial differential model for cardiovascular blood flow is employed. The individual components of the model is represented as an RLC circuit representation. Injury behavior is simulated by varying the circuit system parameters for the section representing different arterial levels. Our results show that significant changes in the blood flow inside the cardiovascular system can be observed when different injury conditions are modeled.

Viscoelastic fluid response can increase the speed and efficiency of a free swimmer.

Microorganisms navigate through complex environments such as biofilms and mucosal tissues and tracts. To understand the effect of a complex medium upon their locomotion, we investigate numerically the effect of fluid viscoelasticity on the dynamics of an undulating swimming sheet. First, we recover recent small-amplitude results for infinite sheets that suggest that viscoelasticity impedes locomotion. We find the opposite result when simulating free swimmers with large tail undulations, with both velocity and mechanical efficiency peaking for Deborah numbers near one. We associate this with regions of highly stressed fluid aft of the undulating tail.

Local flaps: a real-time finite element based solution to the plastic surgery defect puzzle.

One of the most fundamental challenges in plastic surgery is the alteration of the geometry and topology of the skin. The specific decisions made by the surgeon concerning the size and shape of the tissue to be removed and the subsequent closure of the resulting wound may have a dramatic affect on the quality of life for the patient after the procedure is completed. The plastic surgeon must look at the defect created as an organic puzzle, designing the optimal pattern to close the hole aesthetically and efficiently. In the past, such skills were the distillation of years of hands-on practice on live patients, while relevant reference material was limited to two-dimensional illustrations. Practicing this procedure on a personal computer [1] has been largely impractical to date, but recent technological advances may come to challenge this limitation. We present a comprehensive real-time virtual surgical environment, based on finite element modeling and simulation of tissue cutting and manipulation. Our system demonstrates the fundamental building blocks of plastic surgery procedures on a localized tissue flap, and provides a proof of concept for larger simulation systems usable in the authoring of complex procedures on elaborate subject geometry.

Fracturing rigid materials.

We propose a novel approach to fracturing (and denting) brittle materials. To avoid the computational burden imposed by the stringent time step restrictions of explicit methods or with solving nonlinear systems of equations for implicit methods, we treat the material as a fully rigid body in the limit of infinite stiffness. In addition to a triangulated surface mesh and level set volume for collisions, each rigid body is outfitted with a tetrahedral mesh upon which finite element analysis can be carried out to provide a stress map for fracture criteria. We demonstrate that the commonly used stress criteria can lead to arbitrary fracture (especially for stiff materials) and instead propose the notion of a time averaged stress directly into the FEM analysis. When objects fracture, the virtual node algorithm provides new triangle and tetrahedral meshes in a straightforward and robust fashion. Although each new rigid body can be rasterized to obtain a new level set, small shards can be difficult to accurately resolve. Therefore, we propose a novel collision handling technique for treating both rigid bodies and rigid body thin shells represented by only a triangle mesh.

Dynamic simulation of articulated rigid bodies with contact and collision.

We propose a novel approach for dynamically simulating articulated rigid bodies undergoing frequent and unpredictable contact and collision. In order to leverage existing algorithms for nonconvex bodies, multiple collisions, large contact groups, stacking, etc., we use maximal rather than generalized coordinates and take an impulse-based approach that allows us to treat articulation, contact, and collision in a unified manner. Traditional constraint handling methods are subject to drift, and we propose a novel prestabilization method that does not require tunable potentially stiff parameters as does Baumgarte stabilization. This differs from poststabilization in that we compute allowable trajectories before moving the rigid bodies to their new positions, instead of correcting them after the fact when it can be difficult to incorporate the effects of contact and collision. A poststabilization technique is used for momentum and angular momentum. Our approach works with any black box method for specifying valid joint constraints and no special considerations are required for arbitrary closed loops or branching. Moreover, our implementation is linear both in the number of bodies and in the number of auxiliary contact and collision constraints, unlike many other methods that are linear in the number of bodies, but not in the number of auxiliary constraints.

Creating and simulating skeletal muscle from the visible human data set.

Simulation of the musculoskeletal system has important applications in biomechanics, biomedical engineering, surgery simulation, and computer graphics. The accuracy of the muscle, bone, and tendon geometry as well as the accuracy of muscle and tendon dynamic deformation are of paramount importance in all these applications. We present a framework for extracting and simulating high resolution musculoskeletal geometry from the segmented visible human data set. We simulate 30 contact/collision coupled muscles in the upper limb and describe a computationally tractable implementation using an embedded mesh framework. Muscle geometry is embedded in a nonmanifold, connectivity preserving simulation mesh molded out of a lower resolution BCC lattice containing identical, well-shaped elements, leading to a relaxed time step restriction for stability and, thus, reduced computational cost. The muscles are endowed with a transversely isotropic, quasi-incompressible constitutive model that incorporates muscle fiber fields as well as passive and active components. The simulation takes advantage of a new robust finite element technique that handles both degenerate and inverted tetrahedra.